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I review the strategies which have been developped in recent years to solve the non-perturbative renormalization 
problem in lattice field theories. Although the techniques are general, the focus will be on applications to lattice 
QCD. I discuss the momentum subtraction and finite volume schemes, and their application to scale dependent 
renormalizations. The problem of finite renormalizations is illustrated with the example of explicit chiral symmetry 
breaking, and I give a short status report concerning Symanzik's improvement programme to O(a). 



1. INTRODUCTION 

A Quantum Field Theory is renormalized by 
defining the infinite cutoff or continuum limit of 
the regularized theory. The procedure is famil- 
iar in perturbation theory, but is more generally 
applicable. In particular, the renormalization of 
asymptotically free theories at low scales requires 
a genuinely non-perturbative approach. Unfortu- 
nately there exist very few analytical tools, and 
numerical simulations of the lattice regularized 
theory are often the only method to obtain quan- 
titative results. 

One is thus led to discuss renormalization of 
the lattice field theory and to take the continuum 
limit based on numerical data. In general this is 
only possible by making assumptions about the 
continuum approach. Motivated by Symanzik's 
analysis of the cutoff dependence in perturba- 
tion theory, one usually assumes that the con- 
tiuum limit is reached with power corrections in 
the lattice spacing a, possibly modified by loga- 
rithms. In this context "unexpected results" have 
been presented recently by Hasenfratz and Nie- 
dermayer jf]||]. They find that the continuum 
approach of a renormalized coupling in the 0(3) 
non-linear sigma model is compatible with being 
linear in a, rather than quadratic, as one would 
expect for a bosonic model. While the statistical 
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precision of the data is quite impressive, it seems 
premature to draw conclusions. We also note that 
there are many results in (quenched) QCD which 
are well compatible with expectations, and some 
of them will be presented below. Finally, I will as- 
sume the '"standard wisdom" concerning asymp- 
totic freedom and the operator product expan- 
sion (OPE) to be correct. While this has never 
been established beyond perturbation theory, a 
numerical check of the OPE in the 0(3) model 
has recently been presented in ref . || . For an un- 
conventional point of view regarding asymptotic 
freedom I refer to ref. Q . 

Despite the more general title I will focus on 
lattice QCD. For one, non-perturbative renormal- 
ization techniques in QCD have reached a mature 
stage, so that systematic errors can be discussed. 
Second, QCD is a typical example, and the tech- 
niques carry over essentially unchanged to other 
theories. 

This article is organized as follows. In sect. 2, 
I sketch the non-perturbative renormalization 
procedure for QCD. Then the two classes of 
schemes are reviewed which have been proposed 
to solve the problem of scale-dependent renor- 
malization, namely the momentum subtraction 
schemes (RI/MOM, sect. 3) and finite volume 
schemes derived from the Schrodinger functional 
(sect. 4). This is followed by a discussion of finite 
renormalizations (sect. 5), and O(a) improvement 
(sect. 6). I conclude with a remark concerning 
power divergences. 
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2. NON-PERTURBATIVE 

RENORMALIZATION OF QCD 

2.1. Determination of fundamental param- 
eters 

The free parameters of the QCD action are the 
bare coupling and the bare quark mass parame- 
ters. At low energies the theory is usually renor- 
malizcd in a hadronic scheme, i.e. one chooses 
a corresponding number of hadronic observables 
which are kept fixed as the continuum limit is 
taken. 

Once the theory has been renormalized in this 
way, no freedom is left and the renormalized pa- 
rameters in any other scheme are predicted by 
the theory. Of particular interest is the relation 
to the renormalized coupling and quark masses 
in a perturbative scheme. This is equivalent to a 
determination of the A-parameter and the renor- 
malization group invariant quark masses. For, 
given the running coupling g and quark masses 
Mi at scale /i, one finds the exact relations 
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where i labels the quark flavours. Here, (3 and r 
are renormalization group functions with asymp- 
totic expansions, 

d{g) ~ -b g 3 -b ig 5 + ..., (3) 
r{g) ~ </,//" '/:</ ; ' •••• (4) 

Note that the renormalization group invariant 
quark masses are scheme independent, whereas 
the Lambda parameter changes according to 

A' = Aexp{cW/&o}. (5) 

Here, cW is the one-loop coefficient relating the 
renormalized coupling constants g' = g + c^g 3 + 
0(g 5 ). A in the MS scheme and Mi are referred 
to as the fundamental parameters of QCD. 

In order to relate the hadronic scheme to 
these parameters one may introduce an inter- 



mediate renormalization scheme involving quan- 
tities which can be evaluated both in perturba- 
tion theory and non-perturbatively. An obvious 
possibility is to take a renormalized coupling and 
renormalized quark masses which are defined be- 
yond perturbation theory. In such a scheme also 
the renormalization group functions p and r are 
defined non-perturbatively, so that the formulae 
(Q,^) can be applied at any scale /i, once the 
matching to the hadronic scheme has been per- 
formed. 

2.2. Renormalization of composite opera- 
tors 

In the Standard Model composite operators of- 
ten arise when the OPE is used to separate widely 
different scales in a physical amplitude. While the 
operators are originally defined at high energies 
(e.g. at /i = Mw) and in a perturbative scheme, 
one is usually interested in their matrix elements 
between hadronic states. This defines a matching 
problem between a perturbative renormalization 
scheme for the operator (e.g. the MS scheme) and 
a hadronic scheme, where a particular hadronic 
matrix element of the operator assumes a pre- 
scribed value. 

The problem may again be approached by 
defining an intermediate renormalization scheme 
which can be evaluated both non-perturbatively 
and in perturbation theory. At high energies one 
may then either convert to a perturbative scheme, 
or one may construct the renormalization group 
invariant operator, in close analogy to the funda- 
mental parameters introduced above. 

2.3. The problem of scales 

The main problem for the lattice approach is 
due to the large scale differences. In order to 
measure hadron masses the hadrons have to fit 
into the (necessarily finite) space-time volume, 
i.e. L^ 1 <C m OT . On the other hand, the matching 
to the intermediate scheme is done via the bare 
parameters at a matching scale [Ao, and one re- 
quires fiQ <C flmaxj where a max denotes the largest 
lattice spacing considered in the continuum ex- 
trapolation of the matching condition. As the 
available lattice sizes L /a are limited, /^o can not 
be too large, and perturbation theory at this scale 



3 



may not be reliable. This ought to be checked 
by tracing the non-perturbative scale evolution of 
the renormalized parameters in the intermediate 
scheme. If the on-set of perturbative evolution is 
observed around a scale /i, one may then apply 
the formulae ([|||) using the perturbative expres- 
sions (|,f§. 

Combining the various requirements, one gets 

lT x < m n , A < fi < a" 1 , (6) 

and one also wants to vary a in order to take the 
continuum limit. A naive estimate then leads to 
ratios L/a — O(10 3 ) which are clearly too large 
to be realised on a single lattice. 

2.4. Shortcut methods 

Given a hadronic scheme at fixed g one may 
be tempted to avoid introducing the intermediate 
renormalization scheme by considering the bare 
coupling and the bare quark masses themselves 
as running parameters defined at the cutoff scale 
a -1 . One may then match directly to, say, the 
MS scheme at scale /i = a -1 , by using bare per- 
turbation theory || . While this method is quick, 
it is also very difficult to estimate the error in a 
reliable way. In particular it is impossible to dis- 
entangle renormalization from cutoff effects, as 
any variation of a is at the same time a variation 
of the renormalization scale. 

3. THE RI/MOM SCHEME 

The RI/MOM scheme [§ (RI for "regulariza- 
tion independent") is a momentum subtraction 
scheme applied in a non-perturbative context 
pr| . It plays the role of an intermediate renor- 
malization scheme in the sense of the previous 
section. 

3.1. The renormalization procedure 

In the RI/MOM scheme one considers corre- 
lation functions with external quark states in 
momentum space. A necessary first step there- 
fore consists in fixing the gauge (usually Landau 
gauge). One then follows the same procedure as 
in perturbation theory: first the full quark propa- 
gator is computed for a set of momenta and then 
used to amputate the external legs of the corre- 
lation function of interest. The resulting vertex 



function is then considered as a function of a sin- 
gle momentum, and one imposes that its renor- 
malized counterpart be equal to its tree-level ex- 
pression at the subtraction point, viz. 

r R ( P )\ p2= ^T( P r ee . (7) 

For example, in the case of a multiplicatively 
renormalizable quark bilinear operator O eq. (0) 
determines the operator renormalization constant 
Zo{go,a^i) up to the quark wave function renor- 
malization. The latter is then determined by con- 
sidering the vertex function of a conserved cur- 
rent Note that the renormalization constant 
is usually assumed to be quark mass independent. 
This can be achieved by imposing the renormal- 
ization condition (Q) in the chiral limit jl5| . 

The RI/MOM scheme may also be used to ob- 
tain non-perturbatively defined renormalized pa- 
rameters. The PCAC and PCVC relations imply 
that renormalized quark masses can be defined 
using the inverse renormalization constants of the 
axial or the scalar density. To define a renormal- 
ized gauge coupling one may allow for external 
gluon states and impose a RI/MOM condition on 
the 3-gluon vertex function 

3.2. Matching to other schemes 

Suppose one is interested in a hadronic ma- 
trix element of some operator O defined in the 
MS scheme. As a first step one relates the MS 
scheme to the RI/MOM scheme by inserting the 
renormalized operator in the renormalization con- 
dition (Q). This determines the finite renormal- 
ization constant s relating the operators in the 
MS and RI/MOM schemes, order by order in the 
renormalized coupling gl^(p). 

Next, one computes the hadronic matrix ele- 
ment M{a) = (h'\0\h) for a sequence of lattice 
spacings and sets 

(sgjgOO) >WA0 = }™Zo(av)M(a). (8) 

where /i is fixed in terms of a hadronic scale and 
all mass parameters are supposed to be renormal- 
ized in a hadronic scheme. Here one has used the 
regularization independence of the RI/MOM con- 
dition to define the l.h.s. of this equation. Pro- 
vided one has previously determined 5^jg(/i) this 
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completely specifies the matrix element in the MS 
scheme to the calculated order in perturbation 
theory. Of course, it is assumed that continuum 
perturbation theory is applicable at the scale fx. 
Whether this is the case can be checked by repeat- 
ing the matching procedure at a different scale //. 
If the change is compatible with perturbative evo- 
lution from /i to // one may use perturbation the- 
ory to run to larger scales and extract the renor- 
malization group invariant matrix element. 



• r p =i/z p 

♦ r s =i/z s 

- Perturbat. Term 

— P.T. + Power Corr. 



Z m (2GeV, p, 1/a) 




Figure 1. Non-perturbative vs. perturbative scale 
dependence of the renormalization constants of 
the scalar and axial density (ref. ]l2|] , improved 
Wilson fermions at /3 = 6.2). 



3.3. Examples 

Typical examples of operator renormalizations 
are given in fig. [I], where the renormalization con- 
stants for the non-singlet axial and scalar densi- 
ties are plotted against the renormalization scale, 
for 0(a) improved Wilson quarks in the quenched 
approximation at f3 = 6.2 Also plotted are 

the curves obtained in continuum perturbation 
theory. As one observes good agreement with 
the perturbative curves, p-dependent cutoff ef- 
fects are small and perturbation theory seems to 
apply at scales around fx ~ l/a(/3 = 6.2). 

In order to determine the light quark masses, 
the JLQCD collaboration renormalizes the 
scalar density using (quenched) Kogut-Susskind 
fermions at 3 /3-values corresponding roughly to 
a variation of the cutoff by a factor of 2 . The 




Figure 2. Renormalization of the quark mass with 
Kogut-Susskind fermions [ |il"| . Z m should be p in- 
dependent if continuum perturbation theory was 
appropriate. 

renormalization condition is extrapolated to the 
chiral limit, and the renormalized quark mass is 
related to the bare mass m through 
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i(Vo,P, 1/a) m 



(9) 



The total renormalization factor Z m consists of 
the non-perturbatively defined renormalization 
constant at scale p and the perturbative 3-loop 
evolution factor relating the scales p and /xr> If 
continuum perturbation theory was adequate be- 
tween the scales fiQ = 2 GeV and p, then Z m 
should be independent of p. This is not the case 
as fig. H demonstrates. However, if the contin- 
uum extrapolation of eq. (|9|) is performed for 2 
different momenta p = 1.8 GeV and p — 2.6 GeV, 
the results are nicely compatible, demonstrating 
that perturbation theory adequately describes the 
scale evolution in this range of scales. In partic- 
ular, the window condition A -C fx seems to be 
satisfied despite the absence of plateaux in fig. ^ 
at fixed f3. 

Further examples can be found in (? 14 1. In 
particular, the RBC collaboration has applied 
this scheme to the renormalization of operators in 
lattice QCD with Domain- Wall fermions HQ. 
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3.4. Merits and problems of the RI/MOM 
scheme 

The RI/MOM scheme provides a versatile 
framework which can be easily adapted to new 
renormalization problems. Due to the regulariza- 
tion independence, perturbative calculations can 
be performed in a continuum scheme, with many 
results being already available. Furthermore, the 
use of external quark states is not only helpful in 
perturbation theory, but also leads to good sig- 
nals in numerical simulations. 

These nice properties come with a price: a non- 
perturbative gauge fixing procedure cannot avoid 
Gribov copies p8| , and the RI/MOM renormal- 
ization conditions may introduce 0(a) effects into 
the otherwise on-shell 0(a) improved theory. Fur- 
thermore, the continuum perturbation theory is 
done in infinite space-time volume and in the chi- 
ral limit, whereas the numerical simulations are 
performed at finite quark masses and in a finite 
volume with periodic boundary condition. To 
satisfy the RI/MOM conditions chiral extrapo- 
lations are required and one ought to check for 
finite volume effects. Finally, the window con- 
dition A <C /i <C a -1 may not always be satis- 
fied. In particular, problems may be expected 
with operators which couple to the pion. In fact, 
the JLQCD collaboration did not use the axial 
density in the quark mass renormalization, as the 
renormalization constant could not be extrapo- 
lated to the chiral limit JllJ]. This behaviour is 
attributed to a contribution of the Goldstone pole 
which should vanish at high enough scales ||. 
For a more detailed discussion of this problem 
the reader is referred to refs. [ ^9[[lO| ]. 

3.5. Modifications of the RI/MOM scheme 

At this conference, Y. Zhestkov (RBC collab.) 
has presented a first attempt to reach higher 
scales [i before matching to perturbation the- 
ory First results look promising, although 
it is not obvious how finite volume, renormal- 
ization and cut-off effects can be disentangled 
systematically. The situation could be simpli- 
fied if the RI/MOM condition was imposed in 
a finite volume at constant /iL. This would de- 
fine a finite volume scheme, and finite-size scaling 
techniques could then be applied to reach higher 



scales. However, this means that all the pertur- 
bative calculations would have to be re-done in 
the finite volume set-up. 

Finally, attempts are being made to imple- 
ment off-shell 0(a) improvement of the quark 
vertex functions and propagators involved in the 
RI/MOM scheme. A short discussion is deferred 
to sect. 6. 

4. FINITE VOLUME SCHEMES 

4.1. Basic requirements 

Finite volume schemes are obtained by using 
the finite extent of space-time to set the renormal- 
ization scale pp| , pi| , i.e. one identifies \x = L^ 1 in 
eq. (^|). Using recursive finite size scaling meth- 
ods one may bridge large scale difference, so that 
eq. (H) reduces to o < I for any single lattice, 
and applicability of perturbation theory, A <C /i, 
is required only for the large scales (i.e. small L) 
covered by the recursion. 

Finite volume schemes can be defined in many 
ways. However, one would like to maintain 
gauge invariance and respect on-shell improve- 
ment. Furthermore, the renormalization condi- 
tions should be easy to evaluate both by numer- 
ical simulation and in perturbation theory. A 
large family of renormalization schemes with all 
of these properties derives from the Schrodinger 
functional (SF) JUJ^ and is referred to as SF 
scheme. 

4.2. The Schrodinger functional 

The SF is the functional integral for QCD on a 
hyper-cylinder, i.e. all fields are L-periodic in the 
spatial directions, and satisfy Dirichlet boundary 
conditions at Euclidean times and T. For the 
quark fields one sets 

^lxo=0= P> P -^lo=T= J, 

^-lxo=0= P> ^ P +L=T= P'i (10) 

with P± = i(l±7o), and the spatial components 
of the gauge potential take prescribed values C k 
and CI. Renormalizability of QCD in this situa- 
tion has been discussed in (2^j2^] , with the result 
that the SF is finite after the usual renormaliza- 
tion of the parameters and multiplicative renor- 
malization of all fcrmionic boundary fields by the 
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same renormalization constant. 

The SF is a gauge invariant functional of the 
boundary fields, Z = Z[p' , p' ,C; p, p,C], and 
correlation functions may be defined as usual 



(O) = Z 



Oe 



-s 



(11) 

fields J p= pi = p = pi = Q 

where O may contain composite operators and 
boundary fields £, which are obtained by differ- 
entiating with respect to the fermionic boundary 
values, e.g. C( x ) ~~ * <V<5/°( X )- An important point 
to notice is that the gauge field boundary condi- 
tion imply that only global gauge transformations 
are allowed at the boundaries. Bilinear boundary 
sources like C(x)r£(y) are therefore gauge invari- 
ant even for x^y. 
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Figure 3. Continuum extrapolation of the step 
scaling function for the SF coupling with Nf = 
—2, using O(a) improved Wilson "bermions" |27| . 

4.3. Renormalized coupling 

A non-perturbatively defined, gauge invariant 
running coupling can be obtained by choosing a 
family of boundary gauge fields depending on a 
parameter r\ and setting p5| 

d n T 



g\L) = ^ . (12) 

Here, T(r)) = - lnZ[0, 0, C\rj); 0, 0, C(rj)], is the 
effective action of the SF with the perturbative 
expansion, 

r = go 2 r + ri + g 2 Q r 2 + .... (13) 



The numerator in eq. (|l^ ) thus ensures the stan- 
dard normalization of a coupling constant. In 
eq. ((l^) T = L is assumed, the boundary gauge 
fields are scaled with L and the quark masses are 
set to zero. As a result the SF coupling is a quark 
mass independent coupling which runs with the 
scale L. Its perturbative relation to the MS cou- 
pling is known to 2-loop order (cf. [ p6| and refer- 
ences therein). 

For a completely different definition of a finite 
volume running coupling in the SU(2) pure gauge 
theory we refer to ref . [^9| . 

4.4. Non-perturbative running 

The coupling itself is an observable in numeri- 
cal simulations, and the theory may now be renor- 
malized in the chiral limit by imposing 



(14) 



with some numerical value u. The coupling at 
scale 2L is then also determined and can be com- 
puted as a function of u as follows. One chooses a 
lattice size L/a and tunes f3 such as to satisfy (|l4|) . 
Then one measures, at the same bare parameters, 
the SF coupling on a lattice with size 2L/a. This 
yields a lattice approximant to the desired func- 
tion, and by repeating these steps several times 
with increasing resolution L/a one may take the 
limit. 



a(u) = lim T,(u, a/L) _ 2 , r . 



(15) 



An example for such an extrapolation in QCD 
with Nf = — 2 (improved Wilson "bermions" |27| l) 
is shown in fig. |[ Repeating the procedure 
with different values for u the step-scaling func- 
tion a(u) can be constructed in some interval 



recursively 



Setting uq 



1,2, 



one computes 



(16) 



In this way the scale evolution of the coupling is 
traced non-perturbatively and can be compared 
to perturbation theory (cf. figure |J). The AL- 
PHA collaboration has determined the running 
in quenched QCD (2^j3l|] and is currently study- 
ing the theories with Nf = 2 and Nf = — 2 quark 
flavours. The corresponding computer programs 
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for the APE1000 machine are being optimized, 
and the leftmost 3 points in fig. [| have in fact 
been produced by the APE 1000 as a warm-up 
exercise. 
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Figure 4. Non-perturbative running of the SF 
coupling (JVf = 0). For comparison also the per- 
turbative 2-loop and 3-loop curves are shown. 

4.5. Operator renormalization 

To renormalize a composite operator in the SF 
scheme one starts by choosing a non-vanishing 
correlation with a boundary source field. In the 
case of the axial density P a — ip^'y^ip we define 



O a = a 6 ]TC(y)75^C(z) 



(17) 



and similarly 0' a with the primed fields. Given 
the correlation functions 

fe{x Q ) = -\{P a {x)O a ), h = -3^(O a O' a ),(18) 

one may then scale all dimensionfull parameters 
with L, and define the renormalization constant 
through the ratio p(J 

Hi 



Zp ^ a/L) - C 7^wr (19) 

at vanishing quark masses. Here the ratio has 
been chosen such that the boundary field renor- 
malization drops out, and the constant c ensures 
Zp = 1 to lowest order perturbation theory. One 
may then define the step scaling function for the 
renormalization constant, 



ap(u) = Urn 



Zp{go, 2L/a) 



g 2 (L)—u 1 m=0 



a->o Zp(g ,L/a) 
The equations, 

Zp(2L)=ap(u)Zp(L), u = f{L) 



(20) 



(21) 



are then again solved recursively, using the previ- 
ous result for the step-scaling function of the SF 
coupling |3l[| . 

4.6. Universality 

The step-scaling functions are obtained in the 
continuum limit and do not depend on the reg- 
ularization. This is illustrated in fig. [|, which 
shows the continuum extrapolation of two step 
scaling functions, defined in complete analogy 
with (Tp. The simulations have been performed 
both with improved and unimproved Wilson 
fcrmions. Due to incomplete improvement resid- 
ual O(a) effects are expected in both cases, 
and the continuum extrapolated results nicely 
agree. p2] . 

In ]33| an attempt is being made to relate the 
bare chiral condensate obtained with Neuberger 
fermions in p4| to the renormalized axial density 
in the SF scheme. If successful, the step scaling 
function ap could be used to extract the renor- 
malization group invariant chiral condensate. 
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0.05 
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Figure 5. The continuum extrapolation of the 
two step-scaling functions of ref. |32j illustrate 
universality of the continuum limit. The lattice 
data were obtained with both O(a) improved and 
unimproved Wilson fermions. 
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4.7. The static-light axial current 

An interesting new development is the appli- 
cation of the SF techniques to QCD with an in- 
finitely heavy "static" quark j3j|. This is moti- 
vated by the fact that the 6-quark is too heavy to 
be treated as a relativistic particle whilst main- 
taining a physically large volume. Applications 
to physics involving the b-quark therefore rely on 
large extrapolations in the quark mass, or on the 
use of effective theories, such as non-relativistic 
QCD or heavy quark effective theory. Controlling 
QCD in the limit of an infinitely heavy 6-quark 
could lead to a major improvement for quantities 
like Fb , as the extrapolation in the quark mass is 
promoted to an interpolation. 

In ref. |55| the SF for the static light theory has 
been defined in close analogy to the relativistic 
case. In particular it was found that the bound- 
ary quark fields are again multiplicatively renor- 
malized. The application of Symanzik's improve- 
ment programme reveals that the static action is 
already 0(a) improved and determines the struc- 
ture of the improved static light axial current. 
In the effective theory the renormalization of the 
axial current is scale dependent and may be per- 
formed by imposing a similar condition as dis- 
cussed above for the axial density. A potential 
problem is caused by the linear divergence stem- 
ming from self-energy corrections to the static 
quark. Fortunately, the contribution of this coun- 
terterm to the static quark propagator is exactly 
known. It is then possible to choose the ratios 
of correlation functions and parameters such that 
both the boundary field normalization and the 
self-energy counterterms cancel exactly. 

The ALPHA collaboration has obtained the 
step-scaling function in the continuum limit, and 
the non-perturbative running of (a matrix ele- 
ment) of the current is shown in figure [| The 
renormalization problem has thus been solved 
and work concerning the matrix element F B t!Lt is 
in progress. 

5. FINITE RENORMALIZATIONS 

In cases where a non-anomalous continuum 
symmetry is broken by the regularization, finite 
(i.e. scale independent) renormalizations are nec- 
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Figure 6. Non-perturbative running of the static- 
light axial current for Nf — 0. For comparison 
also the perturbative curves are shown ifjBf . 

essary to restore the symmetry up to cutoff ef- 
fects. Well-known examples are the space-time 
symmetry, which is broken up to the discrete sym- 
metry group of the space-time lattice, or global 
axial symmetries, which are completely broken in 
lattice QCD with Wilson quarks. In the case 
of chiral symmetry, the standard solution con- 
sists in imposing the validity of continuum chi- 
ral Ward identities in the theory at finite lat- 
tice spacing |36|. The same idea has been put 
foward in the case of supersymmetry by Curci and 
Veneziano []37f , and first numerical and perturba- 
tive results along these lines have been presented 
at this conference |58|,[59) . A similar treatment of 
the 0(4) Euclidean space-time symmetry appears 
difficult in practice, as it requires the construction 
of the energy momentum tensor, which is a non- 
trivial problem in itself [Q. I will focus on the 
more familiar problem of explicit chiral symmetry 
breaking and its consequences. 

5.1. Chiral symmetry and Wilson fermions 

The explicit breaking of chiral symmetry by the 
Wilson term entails additive quark mass renor- 
malization, a non-trivial renormalization of the 
non-singlet axial current, 

(A R ); = Z A (g )A^, (22) 
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and mixing of operators with the wrong nominal 
chirality. For example, the AS = 2 operator 

A5=2 = 5>"^(i-75)d) 2 , (23) 

decomposes in parity even and parity odd parts, 
O AS = 2 = Ovv+aa - 2 Ova . (24) 

which are renormalized differently. While the 
parity odd operator is multiplicatively renormal- 
ized , the parity even operator mixes with four 
other operators of the same dimension, 

4 

O^+aa = Ovv+aa + 5> (go)0\ d =6] , (25) 

i=l 

and the subtracted operator is then multiplica- 
tively renormalized. 

In order to determine the scale independent 
renormalization constants Z^,Zi,Zvv+AA/ZvA, 
etc., Bochicchio et al. Q have suggested to 
impose continuum axial Ward identities (AWI) 
as normalization conditions. A generic contin- 
uum AWI in the theory with N{ mass degenerate 
quarks has the form 

/ d*x((a^»(s) - 2mP a {x)) 0(y)O ext ) 
Jr 

+ ((£10(y))o ext )=0, (26) 

where R is a finite space time region containing 
y, and O ex t is an arbitrary product of fields local- 
ized outside R. A simple example is the PCAC 
relation which is obtained for 0(y) = S(y — x). If 
imposed on the lattice it determines the additive 
quark mass renormalization constant. Consider- 
ing axial variations of the axial current, the AWI's 
determine the normalization constant Za- More 
generally the axial Ward identities determine the 
relative renormalization of operators belonging to 
the same chiral multiplet, e.g. the scalar and axial 
densities, or Ova and Ovv+aa- 

Applying AWI's has become standard in the 
quenched approximation. At this conference, 
R. Horsley has presented a first attempt to de- 
termine Za for Nf = 2 flavours of 0(a) improved 
Wilson fermions J42| . 

The AWI's can also be implemented more indi- 
rectly by applying the RI/MOM scheme to finite 



renormalization constants. With the appropri- 
ate choice for the renormalization conditions the 
AWI's are automatically satisfied at large enough 
renormalization scales [fl3) . 

5.2. Systematic errors and O(a) ambigui- 
ties 

There is an infinity of axial Ward identities, 
and the results for a given scale independent 
renormalization constant will differ at 0(a) or 
0(a 2 ), if the theory is 0(a) improved. One 
may be tempted to assign a systematic error to 
the renormalization constant corresponding to a 
"typical spread" of the results. This is somewhat 
subjective and unsatisfactory, in particular when 
the spread is large. 

An alternative point of view has been advo- 
cated in refs. @,|]|. One defines a given fi- 
nite renormalization constant Z by imposing the 
Ward identity at fixed renormalized parameters. 
As a result one obtains a smooth function of the 
bare coupling, Z(go), with relatively small errors. 

When Z(go) is used in a continuum extrapola- 
tion the 0(a) ambiguity disappears in the limit. 
This can also be checked explicitly by choosing 
another renormalization condition which yields a 
different function Z(go). The difference Z(go) — 
Z(go) must then vanish in the continuum limit, 
with an asymptotic rate proportional to a or a 2 . 
In ref. p5| the ALPHA collaboration has con- 
firmed this expectation with the example of the 
finite renormalization constant Z m Zp/ZA- 

5.3. Avoiding finite renormalizations 

Some of the scale-independent renormaliza- 
tions with Wilson type quarks can be avoided by 
working in a different basis of fields. Such an 
approach has been presented at last year's con- 
ference |jj6[ , and has been dubbed twisted mass 
QCD (tmQCD). One considers Nf = 2 quark 
flavours with the lattice action, 

= a 4 ^ ip(Dw +m + z/iq7 5 T 3 )-0. (27) 

X 

Here, Dw is the massless Wilson-Dirac operator, 
nio is the standard bare mass, and /x q is referred 
to as twisted mass parameter. This action has al- 
ready appeared in the literature in different con- 
texts • The new idea consists in renormal- 
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izing the theory at finite twisted mass, and in a re- 
interpretation of the renormalized theory as stan- 
dard QCD with Nf = 2 mass degenerate quarks. 
If the renormalization scheme is chosen with care, 
the formulae relating renormalized tmQCD and 
standard QCD can be inferred from classical con- 
tinuum considerations jl9). 

We thus consider the classical continuum la- 
grangian 

C f = i> (p+ m + i/i?75T 3 ) ip + s (p+ m s ) s. (28) 

where ip is a doublet of light quarks, and we have 
also included the strange quark. A chiral (non- 
singlet) rotation of the doublet fields, 

ip' = exp(ia75r 3 /2)-0, 

^' = ^exp(ia7 5 T 3 /2), (29) 

with tana = fi q /m q transforms the Lagrangian 
to its standard form, 

C' f (x)=ij}(x) {p+m')ip'{x) + s(p+m a )s, (30) 

with the light quark mass m' — (m 2 +/iq) 1 ^ 2 . The 
rotation of the quark and anti-quark fields also 
induces a rotation of the composite fields. For 
quark bilinears containing only the light quarks 
one finds e.g. 

A'l = cos(a)4+sin(a)y A 2 , (31) 
= cos(a)^ + 2isin(a)P 3 , (32) 

and the AS = 2 operator behaves as follows, 

°vv+aa = cos(a) Ovv+aa - 2i sin(a) Ova- (33) 

The fields in the primed basis refer to the stan- 
dard QCD basis, and the primed composite oper- 
ators are hence interpreted as usual. Working in 
the twisted basis at a = tt/2 one may use the vec- 
tor current to compute F n , obtain the chiral con- 
densate from P 3 and compute the K — K mixing 
amplitude using the multiplicatively renormaliz- 
able operator Ova jy]- Work on tmQCD by the 
ALPHA collaboration is in progress, and the re- 
sults of a scaling test have been presented at this 
conference [p0[ . 

A variant of the above proposal concerning 
K — K mixing can also be realized with standard 



Wilson fcrmions pHl , by using the axial Ward 
identity (|26|). One considers an axial variation 
of = O VA , 

^A^VA = — |OvV+AA, (34) 

and chooses O ex t as the product of interpolating 
fields for K and K . Imposing the AWI on the lat- 
tice is equivalent to defining the correlation func- 
tion involving Oyv+AA through the remainder of 
the AWI, which contains the operator Ova- The 
trick here is that one directly obtains the desired 
matrix element by choosing C xt appropriately. 
As in the tmQCD proposal this avoids solving the 
complicated mixing problem (|25|), and one is left 
with a scale dependent multiplicative renormal- 
ization for Ova- First feasibility tests have been 
performed and presented at this conference |j52f . 

5.4. Ginsparg- Wilson and Domain-Wall 
quarks 

Massless Ginsparg Wilson fermions are implic- 
itly characterized by a Dirac operator which sat- 
isfies the Ginsparg- Wilson relation j53j 

D lb + ~f b D = aD l5 D. (35) 

The Ginsparg Wilson (GW) relation implies an 
exact chiral symmetry ]54| | , and hence none of the 
finite renormalization constants discussed in this 
section is needed p5[ . 

An explicit solution of the Ginsparg- Wilson re- 
lation has been given by Neuberger |5(|. Unfor- 
tunately, it is computationally very demanding to 
implement GW quarks exactly, and the question 
arises whether an approximation may be good 
enough in practice. A popular choice are Do- 
main Wall fermions p7| , which are formulated in 
five dimensions. The Dirac operator of the cor- 
responding four-dimensional effective theory be- 
comes an exact solution of the Ginsparg- Wilson 
relation as the number of lattice sites in the fifth 
dimension tends to infinity. As the approach is 
expected to be exponential, one may hope that a 
small number N s of points in the fifth dimension 
may already be sufficient. 

To check the size of chiral symmetry viola- 
tions one may pursue the same strategy as in 
the Wilson case, i.e. one measures to what ex- 
tent chiral continuum Ward identities are vio- 
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Figure 7. Difference between the axial and the 
vector current normalization constants as a mea- 
sure of chiral symmetry violation f59fl . 

lated. For example, one may measure the dif- 
ference between vector current and axial current 
renormalization constants The CP-PACS collab- 
oration has measured both renormalization con- 
stants using Ward identities at my = 0, on a 
lattice of size 8 3 x 16 x 16 (cf. fig. @). The 
authors use the methods of rcf. and there- 
fore impose SF boundary conditions in the phys- 
ical time direction. Note, however, it that it 
is not obvious whether this is equivalent to im- 
posing SF boundary conditions in the effective 
4-dimensional theory. The Riken-Brookhaven- 
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Figure 8. The wrong chirality mixing coeffi- 
cients for the AS* = 2 operator with domain 
wall fermions, as determined with the RI/MOM 
scheme Q. 



Columbia (RBC) collaboration has studied the 
mixing properties of the operator Ovv+aa [M . 
They follow ref . pof and determine the wrong chi- 
rality mixing coefficients in a RI/MOM scheme. 
The result for N s — 16 is shown in fig. ||, and one 
observes that the coefficients are indeed small, 
i.e. a factor of 5-10 smaller than the correspond- 
ing result with (unimproved) Wilson quarks J6l| . 

The overall impression is that chiral symme- 
try violations seem to be numerically small in 
many cases. However, it is also possible to deter- 
mine the exponent of the exponential approach 
to a Ginsparg- Wilson solution directly by solv- 
ing a generalized eigenvalue problem |62|]6"5| ] . It 
turns out that the exponent may be rather small 
at typical values of (3, and with Wilson's plaque- 
tte action for the gauge fields. This is potentially 
dangerous, and some modification of the origi- 
nal formulation of Domain Wall Fermions may 
be necessary (cf. |58j for further details). 

6. O(a) IMPROVEMENT 

The purpose of Symanzik improvement ]64| - 
p6| is to accelerate the continuum approach of 
renormalized quantities by introducing appropri- 
ate counterterms to both the action and the com- 
posite operators of interest. The counterterm 
structure follows from the symmetries of the regu- 
larized theory. Without loss of physical informa- 
tion improvement may be restricted to on-shcll 
quantities, i.e. particle masses, energies and cor- 
relation functions of local operators which keep 
a physical distance from each other. In this case 
the equations of motion may be used to reduce 
the number of counterterms. 

6.1. On-shell O(a) improvement with Wil- 
son quarks 

The on-shell O(a) improved action for Nf mass 
degenerate Wilson quarks is given by |^57| , 

Sf = fl4 X! ^( D w + m o + Cswaja^F^ip, (36) 

X 

where Dw denotes the standard Wilson-Dirac op- 
erator. Other counterterms can be absorbed in 
the redefinition of the bare parameters. With the 
subtracted bare quark mass m q = mo — m c , these 
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are |38| 

9l = 



m q (l + o m am q ), 



(37) 
(38) 



and renormalizcd 0(a) improved parameters in 
a mass independent renormalization scheme take 
the form 

9k = 9o z e(.9o,afi), (39) 
tor = fh q Z m (gl,an). (40) 

Improved composite operators are treated simi- 
larly. We restrict attention to non-singlet quark 
bilinear operators, e.g. the axial current, 

(A R )« = Z A (l + b A am q ) {Al + CAad^} . (41) 

The vector current and the tensor density have 
similar additive counterterms, and fields are mul- 
tiplicatively rescaled with a quark mass depen- 
dent term and the appropriate ^-coefficient. Note 
that all renormalization constants are functions 
of the improved bare coupling go, whereas the 
improvement coefficients depend on go. 

6.2. O(a) improvement and chiral symme- 
try 

It is an important observation that all of the 
above counterterms conflict with chiral symme- 
try, and thus repair cutoff effects which origi- 
nate from explicit chiral symmetry breaking. For 
this reason, lattice QCD with Ginsparg- Wilson 
quarks is on-shcll 0(a) improved, provided the 



mass term is introduced in the correct way |69|. 

Concerning Wilson quarks it is a natural ques- 
tion whether axial Ward identities can be im- 
posed as an improvement condition to determine 
the coefficients non-perturbatively. For those co- 
efficients which already appear in the chiral limit 
(the c-coefficients) , this is indeed the case: the 
integral over the region R in eq. (|26| ) can be con- 
verted into a surface integral over OR, and the 
axial Ward identity becomes an identity involving 
on-shell correlation functions only. Away from 
the chiral limit the answer is more subtle. The 
Ward identity now contains an off-shell correla- 
tion function, which is not expected to be im- 
proved. 

An exception is the PCAC relation, which leads 
to an alternative definition of a renormalized 0(a) 



improved quark mass (j40j). Requiring equality 
of the two mass definitions then determines the 
combination 



b m - bp 



9a 



2 9lnZ , 



d9l 



(42) 



where Z = Z m Zp / Z A is a finite renormaliza- 
tion constant. Another special case is the coeffi- 
cient oy which may be determined from the vec- 
tor Ward identities. Alternatively one may define 
an 0(a) improved current starting from the con- 
served current of the exact flavour symmetry, for 
which one finds Zy = 1 and by = 0. 

It is natural to ask whether the situation can 
be improved by allowing for mass non-degenerate 
quarks. Unfortunately the general countertcrm 
structure becomes very complicated, i.e. there are 
many more possible ^-coefficients [[70). Neverthe- 
less, in the case Nf > 3 these authors conclude 
that, assuming that 6 g is known, all but three 
combinations of 6-coefficients are indeed deter- 
mined by Ward identities. 

An independent proposal to determine the 6- 
coefficients is based on the idea that chiral sym- 
metry is restored at short distances Q . This im- 
plies that the physical quark mass dependence is 
invariant under a sign change of the mass. As the 
om q cutoff effects violate this invariance, one may 
in principle be able to isolate the coefficients. Un- 
fortunately, in order to resolve sufficiently short 
distances one may need very small lattice spac- 
ings so that the method may not apply in the 
interesting range of /3-values. 

The situation becomes much simpler if one re- 
stricts attention to the quenched approximation. 
The coefficient b g then vanishes and an improved 
bare quark mass parameter fh q = m q (l + 6 m ar7i q ) 
is defined for each flavour individually. Further- 
more, the flavour off-diagonal quark bilinear fields 
are improved with the same 6-coefficients as in 
the degenerate case, but with m q replaced by the 
average of the valence quark masses in the oper- 
ator |72],(73|]. As shown in ref. |73) it is then in- 
deed possible to determine all 6-coefficients from 
axial Ward identities. In particular, off-shell cor- 
relation functions in the axial Ward identity are 
avoided by considering a partial chiral limit with 
two massless flavours and a further quark mass 
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which is accompanied by the ^-coefficient of in- 
terest. 

6.3. Results 

In perturbation theory the improvement coef- 
ficients are known to one- loop order To 
this order the coefficients do not depend on Nf 
except for b g | fT6| . Nonperturbative results have 
been obtained mostly in the quenched approx- 
imation with the exception of c sw (Nf — 2) [ |78[ . 
Numerical results exist for c sw |7^,|8(| , ca, &v [ p4[ 
and cy |8l|] ■ All these results were obtained by ex- 
ploiting chiral Ward identities formulated in the 
Schrodinger functional, and in general the results 
are obtained as functions of go, for (3 > 6.0. In 
ref. ||72|| the combinations 6a — &p and 6 m were 
determined for values (3 > 6.2, by considering the 
PCAC relation with non-degenerate quarks. Re- 
cently Bhattacharya et al. [g2| have published the 
final version of their work with numerical results 
for all improvement coefficients needed for the im- 
provement of quark bilinear operators. The nu- 
merical results are obtained at two (3 values, 6.0 
and 6.2, using hadronic correlation functions. An 
example is given in fig. ^. 
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Figure 9. The combination oa — by at (3 — 6.2 is 
determined by the slope Q . 

In general, the agreement with previous re- 
sults is reasonably good, except for cy and ca 
at (3 = 6.0. In particular, their preferred value 
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Figure 10. The improvement condition for the 
coefficient ca at (3 = 6.0 [p2[. 



for ca is about —0.04 which is roughly half of the 
value obtained in Q . Fig. 10 illustrates the sen- 
sitivity of the method: It shows the PCAC mass 



obtained from 2-point functions as a function of 
time separation. The excited states in the pion 
channel are sensitive to ca which is tuned such as 
to extend the plateaux to smaller times. This pro- 
cedure is repeated for various quark masses and 
an extrapolation to the chiral limit is performed. 
Similar results for ca with essentially the same 
method have been obtained in ref. Q. 

Improvement coefficients are intrinsically am- 
biguous by terms of O(a). While a large O(a) 
ambiguity in the case of ca cannot be excluded, 
it may also be that the chosen improvement con- 
dition is afflicted by large higher order lattice 
artifacts, resulting in an artificial enhancement 
of the O(a) ambiguity. In ref. j79| some checks 
were carried out to ensure that this does not hap- 
pen and some alternative improvement conditions 
have corroborated the previous result for ca at 
(3 = 6.0 @. 

The ALPHA collaboration has extended the 
work of j72| and determined b m and &a — bp cov- 
ering a range of values (3 > 6.0. Surprisingly, the 

— bp coefficient shows a rather large O(a) am- 
biguity. It was therefore decided to define the im- 
provement condition at "constant physics" , keep- 
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Figure 11. Simulation results and fit function for 
bA — bp, determined from an improvement condi- 
tion imposed at constant physics ]45| . 

ing all quark masses and scales fixed in units of a 
physical scale (cf. sect. 5). The result for the fi- 
nal choice of the improvement condition is shown 
in figure |ll|. As expected the difference to an al- 
ternative definition decreases with a rate roughly 
cx a (cf. fig. Il2|). 



6.4. Off-shell improvement of the gauge 
fixed theory 

In order to render the RI/MOM scheme com- 
patible with 0(a) improvement, one needs to ex- 
tend improvement to the gauge fixed theory (Lan- 
dau gauge), and to off-shell quantities, such as the 
vertex functions and quark propagators entering 
the RI/MOM renormalization conditions. 

According to the authors of ref. Q (cf. 
also Q), this may be achieved by using the 
on-shell improved action supplemented with a 
gauge fixing term and the action of the ghost 
fields. Furthermore, one introduces additional 
0(a) counterterms to the composite operators 
and the quark and anti-quark fields. In the case 
of gauge invariant quark bilinear operators the 
only additional counterterms are the ones which 
vanish by the equations of motion. Improvement 
of the quark field also requires the introduction 
of a non-gauge invariant counterterm oc a^i\> and 
an analogous term for the anti-quark field. 

The latter counterterm is not needed in ex- 
plicit one- loop calculations carried out in ref. |87) . 



Figure 12. The difference between two different 
determinations of 6a — bp at constant physics 
is plotted against a/rQ. One observes the ex- 
pected decrease as the continuum limit is ap- 
proached Ej|. 



These authors organize the calculation differ- 
ently, by explicitly subtracting contact terms in 
the quark propagator and the vertex functions. 
However, this seems to be equivalent to the ef- 
fect of those counterterms which vanish by the 
equations of motion. Finally we mention that 
some of the additional improvement coefficients 
for quark bilinear operators have been determined 
non-perturbatively in ref. |p2|. 



7. CONCLUDING REMARKS 

In the last few years a whole arsenal of tools has 
been created which will help to ultimately solve 
the non-perturbative renormalization problem of 
QCD and similar theories. In this talk I have 
tried to give an impression of the problems and 
the conceptual and technical progress in the field. 

An important omission in this review concerns 
power divergences. The reason is that power di- 
vergences do not only present a technical but 
above all a serious conceptual challenge. To illus- 
trate the potential problem, assume that QCD is 
regularized on the lattice with Ginsparg- Wilson 
quarks, and suppose we want to renormalize the 
isosinglet scalar density at non-zero quark mass. 
The structure of the renormalized operator is de- 
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termined by the symmetries, 

S° =Z S {S° + ma- 2 c s (am,g )} , (43) 

where cs is a dimensionless coefficient which de- 
pends on the bare parameters and is independent 
of the renormalization scale \x jSg] . To define the 
renormalized operator one should be able to im- 
pose a renormalization condition which does not 
refer to the regularization. In the chiral limit 
this is not a problem, as the additive renormal- 
ization vanishes exactly. At non-zero mass, how- 
ever, one must first define a renormalized quark 
mass, which is only possible up to an intrinsic 
0(a 2 ) ambiguity. The renormalization condition 
for the scalar density now refers to the renormal- 
ized quark mass, and one may be worried that 
the 0(a 2 ) ambiguity combines with the quadratic 
divergence to produce an 0(1) ambiguity in the 
"renormalized operator" . 

The example of the isosinglet scalar density 
may seem academic, but similar problems are 
expected in the renormalization of the effective 
weak hamiltonian. A theoretical solution to 
this problem might be Symanzik improvement to 
higher orders. However, this requires the intro- 
duction of four-quark operators in the action and 
does not seem very practical. New ideas may be 
necessary to either solve or circumvent these diffi- 
cult renormalization problems, and an interesting 
new approach has already been proposed [^9| . 

I would like to thank the conference organizers 
for the opportunity to give this talk. I have bene- 
fitted from discussions with many participants at 
the conference, and in particular with M. Luscher, 
R. Petronzio, G.C. Rossi and R. Sommer. Fur- 
thermore I thank R. Sommer for critical com- 
ments on a draft of this writeup. Support by 
the European Commission under grant No. FM- 
BICT972442 and by CERN is gratefully acknowl- 
edged. 
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